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INTERLAMINAR FRACTURE TOUGHNESS: THREE-DIMENSIONAL FINITE ELEMENT 


MODELING FOR END-NOTCH AND MIXED-MODE FLEXURE 

P.L.N. Murthy* and C.C. Cham1s+ 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, ' o 44135 


ABSTRACT 

A computational procedure Is described for evaluating End-Notch-Flexure 
(ENF) and Mixed-Mode-Flexure (MMF ) Interlaminar fracture toughness In unidirec- 
tional fiber composites. The procedure consists of a three-dimensional finite 
element analysis In conjunction with the strain energy release rate concept and 
with composite micromechanics. The procedure Is used to analyze select cases 
of ENF and MMF. The strain energy release rate predicted by this procedure Is 
In good agreement with limited experimental data. The procedure Is used to 
Identify significant parameters associated with Interlaminar fracture tough- 
ness. It Is also used to determine the critical strain energy release rate and 
Its attendant crack length In ENF and/or MMF. This computational procedure has 
considerable versatlllty/generallty and provides extensive Information about 
Interlaminar fracture toughness In fiber composites. 


INTRODUCTION 

Interlaminar delamlnatlon of composites Is a type of fracture mode which 
needs to be carefully examined and properly considered In the design of com- 
posite structures. Regions prone to delamlnatlons Include free edges, loca- 
tions of stress concentration. Joints, Inadvertent damaged areas and defects 
resulting from the fabrication procedure. 

One way of properly accounting for Interlaminar delamlnatlon In a design 
Is to determine Interlaminar fracture toughness parameters and then to evaluate 
those stress states which are likely to Induce Interlaminar fracture. Several 
test methods to determine fracture toughness have been proposed and are cur- 
rently being used. These tests Include: edge delamlnatlon, double-cantilever 

beam (constant and variable thickness), cracked-lap-shear, biaxial Interlaminar 
fracture (ARCAN), and the recently Introduced three-point bend tests for Mode 
II (end-notch flexure (ENF)) and mixed Modes I and II fracture (MMF). Each of 
these test methods has Its advantages and limitations (figs. 1 and 2). These 
tests have been the subject of discussion and evaluation at ASTM D30.02 and 
D30.04 subcommittee meetings and specialty symposia sponsored by these sub- 
committees (refs. 1 to 3). 

At this time It appears that the three-point bend tests, (ENF and MMF, 
fig. 2) have some unique features over the others-especlal ly for determining 
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Interlaminar Mode II and mixed mode fracture toughness. These unique features 
stem from (1) the simplicity of the test, and (2) the ability to measure the 
fracture toughness parameters directly during the test. Recent research 
efforts at Lewis Research Center have focused on the development of a computa- 
tional method (procedure) for simulating the three-point bend test and evalu- 
ating Interlaminar fracture toughness In unidirectional fiber composites as 
determined by ENF and/or MMF. This computational procedure consists of a 
three-dimensional finite element analysis In conjunction with the strain energy 
release rate concept and with composite mlcromechanlcs . The procedure Is suit- 
able for determining global and local Interlaminar fracture toughness param- 
eters as well as critical values of these parameters. The objective of this 
report Is to describe the computational procedure In detail and results 
obtained therefrom. 

A unique feature of this computational procedure Is the consideration of 
the Interply layer as a distinct entity and with finite thickness. The Inter- 
ply layer has not previously been Included In three-dimensional laminate 
analysis (ref. 4) and In fracture mechanics-type computations. The Interply 
layer needs to be considered In order to compute accurately those stress fields 
which Induce Interply delamlnatlon. 


FRACTURE TOUGHNESS PREDICTIONS 

MSC/NAS7RAN three-dimensional finite element static analysis (FEA) with 
substructuring was used to determine the structural response variables (dis- 
placements, stresses, strains) required for fracture toughness predictions. 
These structural response variables were used subsequently with three different 
methods to predict the Interlaminar and mixed mode fracture toughness. The 
three methods used are: (1) the global method, (2) the local crack closure 

method, and (3) the NASA Lewis "unique" local crack closure method developed 
during this Investigation. Each method Is summarized below. 


Global Method 

The specific computational steps for this method are as follows (refer to 
f 1 o d). 

(1) Model the specimen with crack length (a) using three-dimensional fin- 
ite elements as described later. 

(2) Apply a load (P; at specimen midspan. 

(3) Calculate the midspan displacement v(a) using three-dimensional FEA 
as described later. 

(4) Induce crack extension Aa keeping load (P) constant. 

(5) Calculate the midspan deflection [v(a + Aa)]. 

(6) Determine the Strain Energy Release Rate (SERR), G, from 
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G - P x [v(a ♦ &a) - v(a)j/2b&a 


( 1 ) 


where b Is the specimen width. 

(7) Repeat steps (4) to (6). 

(8) Plot results for G versus a or Aa. 

(9) Identify fracture toughness characteristics as described later. 

(10) Examine complete stress state near crack tip. 

(11) Conpare with corresponding uniaxial composite strengths. 

(12) Look for possible correlation of fracture toughness with conposlte 
uniaxial strengths. 

The global method yields the global fracture toughness without any regard to 
participating and/or dominating local fracture modes. 


Local Crack Closure Method 
The specific steps for this method are as follows: 

(1) Perform steps (1) and (2) ss In the Global Method. 

(2) Calculate (u, v, w) a at the -'■ack tip nodes. 

(3) Induce crack extension Aa keeping P constant. 

(4) Calculate (u, v, w) a + aa at the same nodes as In step (2). 

(5) Apply unit forces (f x , f y , f z ) at these nodes while keeping (P) 
and (a *■ Aa) constant. 

(6) Calculate corresponding (u, v, w) displacements. 

(7) Calculate local forces required to "close" the crack In Its respective 
planes (fig. 2) using the following equations 




(8) Determine the local SERR's from 

Gj = F y x [v(a + Aa) - v(a)]/2b Aa (5) 
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Gj i = x [u(a ♦ Aa) - u(a)]/2bAa 


( 6 ) 
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Gj ii = x [w(a * Aa) - w(a)]/2bAa (7) 

(9) Repeat steps (3) to (8). 

(10) Follow steps (8) to (12) in the Global Method. 

The local crack closure method yields the contribution of each local 
fracture mode to the composite Interlaminar or mixed mode fracture toughness. 
The global fracture toughness can also be determined since the midspan dis- 
placement Is available from the FEA. 


NASA Lewis ’Unique" Local Crack Closure Method 

The specific steps c or this method are as follows: 

(1) Perform steps (1) tc * s; In the Local Crack Closure Method. 

(5) Apply enforced displacements (single point constraints) using the step 
i/i ol splacements (u, v, w^ a at the crack tip nodes. 

(6) Repeat FEA with these single point constraint.. 

(7) Calculate the corresponding forces at these constraints (F x , Fy, F z ). 
These are called the single point constraint forces In FEA. 

(8) Calculate the respective SERRs using equations (5) to (7) with the 
F x , Fy, and F z from step (7). 

(9) Repeat steps (1) to (8). 

(10) Follow steps (8) to (12) in the Global Method. 

The NASA Lewis method is a arlation of the local crack closure method. 

The variation arises from the unique FEA feature: the single point constraint 

force. It is also conceptually simpler than the local crack closure method and 
has the added advantage of preserving the characteristic of monotonically 
increasing displacement under the applied load with crack extension. On the 
other hand, It is possible to obtain reversal in this displacement when using 
the local crack closure method. Comparisons of results predicted using these 
two methods are made In a later section. 


FINITE ELEMENT ANALYSIS MODELS 

The entire specimen was modeled, Including the Interply layer, using the 
three-dimensional finite elements available In MSC/NASTRAN. A schematic of the 
model with the actual dimensions used is shown In figure 3. The specimen 
modeled Is the same as that currently being evaluated by the ASTM D30.04 sub- 
committee for a possible standard test method to determine Mode Ii (using ENF) 
and mixed Mode I and II (using MMF) fracture toughness. The Interply layers 
and the Individual plies are modeled in the vicinity of the crack. The 
remaining plies are grouped as shown schematically in figure 4. 
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A three-dl wens Iona 1 conputer plot of the finite element Model Is shown In 
figure 5. The speclaen was Modeled using 1536 solid eleMents for a total of 
6090 degrees of freedom (OOF). A superelement was used In the crack vicinity. 

A two-d 1 Mens Iona 1 computer plot of the superelement Is shown In figure 6. The 
superelement contained 360 solid elements, 450 nodes (1350 DOF). The displace- 
ment boundary conditions (fig. 3) are: 
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A line load was applied at the specimen midspan (fig. 2). The magnitude 
of this load was 120 lb. This magnitude corresponds to the experimental load 
which Induced unstable crack propagation as reported In the ASTM D30.02 sub- 
committee meeting (April 1984). The crack extension propagation was simulated 
by progressively deleting Interply layer elements and repeating the FEA as 
described In the previous section. The justifications for deleting the element 
Instead of the conventional line opening using double nodes are as follows: 

(1) The Interply layer element Is very thin (about 10 percent of the 
fiber diameter) . 

(2) Photomicrographs of the fractured surface generally show relatively 
little resin residue Indicating that the Interply layer is destroyed during 
the fracture process. 

(3) The crack does not remain at the midplane of the Interply layer but 
meanders between the fiber surfaces of the adjacent plies. 

The composite material simulated was AS-graphlte fiber/epoxy matrix (AS/E) 
unidirectional composite. The constituent material properties are summarized 
In table I assuming room temperature dry conditions. Five different unidirec- 
tional composites were simulated with the following fiber volume ratios (FVR): 
0.30, 0.55, 0.60, 0.65, and 0.75. The unidirectional ply properties, the 
Interply layer thickness and Its properties, and the appropriate material 
properties required for the three-dimensional finite element analysis were 
generated with the aid of ICAN (Integrated Composite Analyzer), (ref. 5). The 
properties for the Interply layers were assumed to be the same as those for the 
matrix. The three-dimensional composite properties (using MSC/NASTRAN designa- 
tion) and Interply layer thickness predicted by ICAN are summarized In 
table II. 
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Even though room temperature dry conditions were assumed In this simula- 
tion, hygrothermal environmental effects can readily be simulated using the 
theory (ref. 6) embedded In ICAN to predict the corresponding composite prop- 
erties. Other composite systems Including (1) Intraply hybrids, (2) Interply 
hybrids, (3) Intra/Inter ply hybrids, (4) composites with perforated Interply 
layers, and (5) composites with different void content can readily be simulated 
using ICAN. 


RESULTS AND DISCUSSION 

The results obtained together with appropriate discussions. Interpreta- 
tions, and possible Implications are summarized below first for the global 
method and second for the local crack closure methods. The results are 
presented graphically where midspan displacement and SERR are plotted versus 
crack propagating length for the five different fiber volume ratios (FVR). In 
addition, graphical results are presented for the stress field behavior In the 
vicinity of the crack and also the variation of maximum stress magnitudes 
versus crack propagating length. The maximum stress magnitudes are assumed to 
be at the centroid of the Interply element In the superelement ahead of the 
crack tip. 


Global Method Results 

The midspan displacement versus crack propagating length for the five 
different fiber volume ratios Is shown In figure 7 for the ENF specimen. . The 
midspan displacement decreases with Increasing fiber volume ratio. Also, the 
midspan deflection Increases gradually with crack opening length up to about 
a =1.20 In, and then Increases rapidly. The corresponding SERR Gjj 
determined using the global method described previously (Eq. (1)) Is shown In 
figure 8. The authors Interpret the behavior shown In figure 8 to be 
associated with Interlaminar composite fracture characteristics as follows: 

(1) The Initial portion of the curves (a < 1.05 In) represents slow crack 
growth. 

(2) The Intermediate portion (1.05 In < a < 1 .20 In) represents ^cable 
crack growth. 

(3) The rapid Increase (rise) (a > 1.20 In) represents unstable, and 
therefore rapid crack propagation to fracture. 

The SERR Gjj varies approximately Inversely with FVR as would be 
expected since It Is calculated from the midspan displacement which also varies 
inversely with FVR. 

The corresponding results for MMF are shown In figure 9 for displacement 
and In figure 10 for SERR G. These curves exhibit the same behavior as those 
for ENF and, therefore, exhibit similar Interlaminar composite fracture 
characteristics. A comparison of the curves In figures 8 and 10 reveals that 
the Interlaminar fracture characteristics for ENF and MMF are almost Identical 
for the same load. Again, this Is anticipated since they are both determined 
using the midspan displacement which Is about the same for the two cases. 
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The results In figures 8 and 10 show that the anticipated "global" com- 
posite Interlaminar fracture strain energy release rate can be predicted by 
conducting three-dimensional finite element analyses on ENF and MMF specimens 
and using the global method. In addition, these results show similar char- 
acteristics for shear (Mode II) and mixed (Mode I and II) fracture. One con- 
clusion from the above discussion Is that the global method Is quite general 
for determining the mixed mode composite Interlaminar fracture characteristics. 
Embedded cracks, multiple cracks, crack location, specimens of uniform and 
variable cross sections, and different laminate configurations can be just as 
easily handled as the two cases already described. The global method generally 
does not separate the contribution of each mode. However, these contributions 
can be determined by the local crack closure methods as described In the next 
section. 


Local Crack Closure Methods Results 

The Mode II SERR, determined using the local crack closure method 
(eq. (6)), Is shown In figure 11 (dashed line) for 0.6 FVR. The curve deter- 
mined using the global method Is also plotted for comparison. The range of 
experimental data reported In the ASTM 030.02 subcommittee meeting Is shown by 
dashed horizontal lines. Three points are worth noting In figure 11: 

(1) The global method predicts higher Gj j values than the local crack 
closure method. 

(2) The two methods predict similar composite Intralaminar fracture 
behavior . 

(3) The measured data coincide with the rapid rise In both curves. 

Accordingly, two general conclusions can be drawn: 

(1) The rapid rise In the SERR (Gn) versus crack length (a) coincides 
with the measured critical Gj j . 

(2) A conservative assessment of the composite Interlaminar fracture 
characteristics can be obtained by using the Gjj determined from the local 
crack closure method. 

The mixed mode (I and II) composite Interlaminar fracture characteristics, 
determined from the local crack closure method (eqs. (5) and (6)), are shown 
In figure 12 for a composite with 0.6 FVR. The fracture characteristics for 
each mode (Mode I and Mode II) are compared with those predicted for the mixed 
mode by using the global method (solid line) and the algebraic sum of the two 
modes from the local method (short dashed line). The following observations 
are worth noting In figure 12: 

(1) Both the local crack closure and the global methods predict similar 
mixed mode (Mode I and II) fracture characteristics. 

(2) The local crack closure method predicts lower mixed mode G values 
than the global method. 
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(3) Mode II fracture dominates the crack extension during slow and stable 
crack growth (a < 1.05 In). 

(4) Mode I dominates the crack extension during unstable crack growth and 
rapid crack propagation (a > 1.15 In). 

(5) Mode I has slightly negative values (a < 1.15 In) Indicating local v 
(predicted by eq. 5) displacement reversal. 

The general conclusions to be drawn from the above observations are: 

(1) The contribution of each fracture mode can be determined using the 
local crack closure method. 

(2) Mode I drives the composite Interlaminar delamlnatlon In HMF 
specimens . 

(3) The local method provides a conservative assessment of mixed mode 
composite Interlaminar fracture characteristics as was the case for Mode II. 

The mixed mode (Mode I and II) composite Interlaminar fracture charac- 
teristics predicted using the Lewis method (local single point constraints) are 
shown In figure 13. T he characteristics predicted by using this method are 
practically the same as those observed In figure 12. One exception Is that 
Mode I Is positive for all 'a' and, thus, this method preserves the v dis- 
placement monotonicity. A typical value for Gj Is about 0.6 ps 1 - in at the 
onset of rapid crack propagation for this type of composite. This coincides 
with the rapid rise In Mode I (figs. 12 and 13). It Is Important to note that 
the contribution of each mode to mixed mode fracture can be jetermlned, within 
acceptable engineering accuracy, by either of the local mr.hods described 
herein. The local crack closure method, however. Is easier to Implement 
computationally based on the authors' research experlenci . 

The composite Interlaminar opening mode fracture was also determined by 
simulating a double cantilever beam (OCB) specimen with the same geometry and 
loads used In the ENF and MMF. The results are shown In figure 14 for both the 
global and the local methods. The general fracture characteristics are similar 
to those determined for Mode I by the local methods of the MMF (figs. 12 and 
13). The significant exception Is that the opening fracture mode of the DC 
specimen becomes unstable at a smaller crack opening compared to the MMF 
specimen . 

Apparently the presence of Mode II Increases the resistance to Mode I 
fracture. This Is consistent with experimental observations (ref. 7). The 
significant conclusion Is that the local crack closure methods predict the 
magnitude of the opening mode contribution and Its dominance to rapid crack 
propagation In mixed mode fracture reasonably well. 


LOCAL STRESS FIELDS 

End-Notch-Flexure . - The averaged across-the-wldth stress field behavior 
In the Interply layer In the ENF Is shown In figure 15. Interestingly, all 
three stresses o xx , o y y , and o xy have approximately the same magnitude In 
the element next to the crack tip. The two normal stresses (o xx and o yy ) 
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are compressive, monotonic, and decay rapidly with distance from the crack tip 
as would be expected from the near-singular stress field. On the other hand, 
the shear stress exhibits oscillatory behavior near the crack tip, decreases 
to a relatively constant value and then slowly decays to the value predicted 
by simple beam theory. At this time the a ithors suspect that this oscillatory 
behavior may be real. Two supporting reasons are: 

(1) The proper behavior of the two normal stresses gives confidence In the 
analysis, and 

(2) The relatively constant value of the shear stress at a distance from 
the crack-tip greater than 0.05 In. 

The peak stress field behavior (stress In the Interply layer crack tip 
element) versus the extending crack length Is shown In figure 16 for 0.6 FVR. 
Also shown In this figure Is the Interlaminar shear strength predicted by using 
the composite micromechanics equations programmed In ICAN. The peak shear 
stress (o xy ) reaches relatively high values during the slow crack growth 
stage and remains practically constant during the stable crack growth stage 
(fig. 8). The peak shear stress exceeds the corresponding strength for 
practically the entire extending crack distance. The peak value of the 
through-the-thlckness normal stress (oyy) Is substantial (comparable to 
o xx ) and It Is compressive throughout tne extending crack distance. This 
compressive normal stress will tend to close the crack and thus r-.hance the 
Interlaminar resistance to Mode II type fracture. 

It appears that the rise In the peak shear stress Is associated with an 
Increase In the normal compressive o y y stress. The peak normal stress 
(o xx ), along the specimen span, Is also compressive and has about the same 
magnitude as oyy. Both of the normal compressive stresses tend to crush the 
fractured Interply layer and, therefore, provide additional credence to sim- 
ulating the extending crack by deleting the Interply layer crack-tip element. 

The significant conclusion from the above discussion Is that Interply Mode 
II fracture occurs when the corresponding shear strength Is exceeded. This 
conclusion can be used to determine critical fracture toughness parameters as 
will be described later. 

Mixed -Mode- Flexure . - The corresponding average stress-field behavior In 
the MMF Is shown In figure 17. The two normal stresses (o yy and o xx ) 
exhibit singular monotonic stress behavior while the shear stress (o X y) 
exhibits oscillatory behavior near the crack tip and becomes constant there- 
after. Three points worth noting are: 

(1) oyy has the highest magnltudt In the crack-tip element. 

(2) Both normal stresses are tensile while they were compressive In the 
ENF specimen. 

(3) The shear stress magnitude beyond the crack tip vicinity Is about the 
same as that In the ENF flexure specimen (fig. 15). 

The peak stresses In the MMF specimen versus extending crack-tip distance 
are shown In figure 18. Also shown In this figure are the corresponding 
Interlaminar shear and transverse tensile strengths predicted using ICAN. 
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Both stresses (oyy and o xy ) exceed their corresponding strengths for 
practically the entire extending crack-tip distance. Additionally, oyy 
rises more i-r.dly while t* X y remains approximately constant for extending 
crack-tip distances greater than 0.150 In. These result are consistent with 
the SERR G behavior (fig. 10) and demonstrate that the MMF specimen Is sub- 
jected to both opening mode (Mode I) and shear mode (Mode II) type fractures. 

The significant conclusion from the above discussion Is that mixed mode 
fracture occurs In MMF specimens when the stresses oyy and/or a xy 
exceed their corresponding strengths. This Is consistent with the conclusion 
stated previously for the ENF specimen. It has the same significant Implica- 
tions and can be used to determine mixed mode fracture toughness as will be 
described later. 


FRACTURE TOUGHNESS PARAMETERS 

During the course of this Investigation, It became evident that composite 
Interlaminar fracture toughness Is associated with several Intrinsic properties 
of the composite. The most obvious ones Include Interply 3 ar thickness. 
Interlaminar shear strength and transverse tensile strength. These are In 
addition to the well known critical strain energy release rate (SERR), critical 
stress Intensity factor and crack length. Furthermore, the Intrinsic composite 
properties depend on composite constituent material properties and on fiber 
volume ratio (FVR) and, as a result, are not Independent. Their Interdepend- 
ence can be determined by composite micromechanics (refs. 8 and 9). 

The effects of fiber volume ratio on the fracture toughness parameters are 
presented graphically In figure 19 for ENF (Mode II) fracture and for a 120 lb 
load. The parameters plotted are: Interlaminar shear strength, SERR and 
Interply layer thickness for three different crack lengths (a a). The Interest- 
ing points to be noted are: 

(1) The SERR decreases with Increasing FVR (rapidly at low FVR (about 
0.30) and more slowly at high FVR (about 0.70)). 

(2) Both the Interlaminar shear strength and the Interply layer thickness 
decrease with Increasing FVR and at a progressively faster rate. 

Both of these lead to the conclusion that fiber composites with relatively 
high FVR will have lower fracture toughness compared to those with Intermediate 
(about 0.5) or low FVR. Another, and perhaps by far more Important, conclusion 
Is that composite mechanics In conjunction with three-dimensional finite ele- 
ment analysis provide unique computational analysis methods for descrlblng/- 
assesslng the fracture toughness of fiber composites In general. 


PREDICTION OF COMPOSITE INTERLAMINAR FRACTURE TOUGHNESS 

The results of this Investigation taken collectively lead to a general 
computational procedure for predicting unidirectional composite Interlaminar 
fracture toughness using the End-Notch-Flexure (ENF) and/or Mixed-Mode-Flexure 
(MMF) methods. In developing this procedure, the critical parameters associ- 
ated with the fracture toughness of the composite are revealed. This procedure 
consists of the following steps: 
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(1) Determine the requisite properties at the desired conditions using 
composite micromechanics. These requisite properties will be similar to those 
summarized In table II. 

(2) Run a three-dimensional finite element analysis on an CNF or MMF 
specimen (fig. 2) for an arbitrary load using finite element models similar to 
those described herein (figs. 3 and 4). Note that the superelement Is not 
necessary although It expedites the computations. 

(3) Scale the arbitrary load In step (2) to match Interlaminar shear 
stress or thrcugh-the-thlckness normal stress In the crack-tip element with the 
Interlaminar shear strength or the through-the-thlckness normal strength of the 
composite. Note that this element Is next to the crack tip. 

( 4 ) Pun several finite element analyses with the scaled load and progres- 
sively extending crack length (a). 

(5) Calculate strain energy release rates (G) using the appropriate 
equations (eqs. (1), (5), (6), or (7)). 

(6) Plot the curve G calculated at step (5) versus a. This curve Is 
similar to those In figures 11 to 13. 

(7) Construct two tangents to the curve as follows: 

(a) through the slowly rising portion of the stable crack growth 
region 

(b) through the rapidly rising portion of the unstable crack growth 
region. 

(8) Select the critical G (strain energy release rate) value and the 
critical crack-length from the Intersection of the two tangents In step (7). 

The steps for this procedure are summarized together with a schematic In 
figure 20 for convenience. 

It Is Important to note that this procedure has not been verified with 
appropriate experimental data other than that described herein. Therefore, It 
should be used judiciously until sufflcently well validated. However, it 
provides a convenient computational method for Initial screening of candidate 
composites to select those which should meet Interlaminar fracture toughness 
design requirements. It also provides a direct means for quantifying values 
for critical strain energy release rate (G) and Its corresponding critical 
length (a) on which design criteria can be based. It Is the authors' con- 
sidered techt.lcal Judgment that proper use of this computational procedure 
provides detailed Information which can be used to describe, assess, and 
predict composite Interlaminar fracture toughness. As a resiult, composite 
structures can be designed to meet this hitherto elusive critical design 
requirement. This computational procedure has considerable versatility and 
generality as was previously discussed. 


SUMMARY OF RESULTS 

The significant results of an Investigation of composite Interlaminar 
fracture determined by the computational simulation of end-notch-'lexure (ENF) 
and mixed-mode-flexure (MMF) are summarized below: 
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1. Three-dimensional finite element analysis In conjunction with composite 
mechanics provides a direct computational method for determining the Inter- 
laminar fracture toughness In unidirectional composites. 

2. The significant composite properties associated with composite Inter- 
laminar fracture toughness are: (a) Interlaminar shear strength, (b) trans- 

verse tensile strength, (c) Interply layer thickness, and (d) fiber volume 
ratio. 

3. Three methods are suitable to computationally determine the Inter- 
laminar strain energy release rate (SERR): (a) the global method, (b) the 

local crack closure method, (c) the local single point constrained force 
(Lewis) method. However, only the local methods are suitable for determining 
the contributions of each fracture mode. 

4. The global method predicts higher SERR compared to local methods during 
slow and stable crack growth. However, all three methods tend to converge 
during rapidly Increasing strain energy release rate (SERR) which Indicates 
unstable crack growth. 

5. Predicted SERR for shear (Mode II) and mixed mode (I and II) Inter- 
laminar fracture are In good agreement with limited measured data. 

6. Mode II (shear) dominates ENF while Mode I (opening) dominates the 
unstable crack growth In MMF. 

7. Mode II (shear) unstable crack growth occurs In ENF when the Inter- 
laminar shear stress intensity (as determined herein) exceeds Its corresponding 
shear strength. 

8. Mixed mode (I and II) unstable crack growth occurs In MMF when the 
through-the-thlckness normal stress Intensity (as determined herein) exceeds 
Its corresponding transverse tensile strength. 

9. A computational procedure was developed to determine the critical SERR 
and Its attendant crack length associated with ENF and/or MMF fracture. 
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TABLE I. - CONSTITUENT MATERIAL PROPERTIES USEO IN THE 

SIMULATION 


Property ! 

| j 

i 

; j 

Units 

Fiber 

AS-graphlte 

Matrix 

Intermediate-modulus 

high-strength 

Symbol 

L _ 

Value 

Symbol 

Value 

• r 

i Elastic moduli j 


wm 

31.00 

Em 

0.5 

j i 


Ef 22 

2.00 





G f T 2 

2. 





G f 23 

1. 

Gm 

.1852 

Poisson's ratio 


Vf 12 

0.2 

Vm 

.35 



v f 23 

0.25 



Strengths 

ksl 

SfT 

400 

SmT 

15 



Sfc 

400 

^mc 

35 





^ms 

13 

Fiber diameter 

Inch 

d f 

— 

0.0003 








TABLE II. - COMPOSITE PROPERTIES PREDICTED BY ICAN AND USED IN THREE-DIMENSIONAL FINITE ELEMENTS 












I 


EDGE-DELAMINATION 
MIXED MODE 
(I& III 



DOUBLE CANTILEVER 
OPENING MODE (II 



figure 1. - Schematics of test methods for measuring interlaminar fracture toughness. 



Figure 2. - Schematic of flexural test for interlaminar fracture mode toughness. 
Note: Origin at left support bottom. 




1536 SOLID ELEMENTS 
(32 PENTAHEDRANS AND 
1504 HEXAHEDRANS) 

6090 DOF 

SUPER ELEMENT IN CRACK 
REGION 


Figure 5. - 3-0 finite element model. 



Figure 6. - Crack region superelement model details-front view. Superelement statistics: 
360 solid elements (32 6-node pentahedrans, 32S 8-node bricks) 450 nodes, 1350 DOF. 






















Figure 17. - Stress field in interply layer near crack tip. Mixed-mode-flexure 
(AS/Ei, 



Figure 18. - Peak stres'. behavior at rrack tip with crack extension. Mixed-mode- 
fracture IAS/EI. 
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Figure 19. - Fiber volume ratio effect on interlaminar fracture 
toughness parameters. End-notch-flexture (AS/EU 


• DETERMINE REQUISITE PROPERITES AT DESIRED CONDITIONS USING COMPOSITE 
MICROMECHANICS 

• RUN 3-D FINITE ELEMENT ANAIYSIS ON ENF IMMFI SPECIMEN FOR AN ARBITRARY LOAD 

• SCALE LOAD TO MATCH INTERLAMINAR SHEAR STRESS AT ELEMENT NEXT TO 
CRACK-TIP 

• WITH SCALED LOAD EXTEND CRACK AND PLOT STRAIN ENERGY RELEASE PATE 
VS CRACK LENGTH 

• SELECT CRITICAL "G " ) 

AND CRITICAL "a" ) 

• METHOD HAS VERSATILITY (GENERALITY Extended crack length, a 

Figure 20. - General procedure for predicting composite Interlaminar frac*ure toughness 
using the end-nctch-flexure IENFI or mixed-mode-flexure (MMF) method. 
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